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tH- ■ The generation of the baryon asymmetry of the universe is con- 

sidered in the standard model of the electroweak theory with simple 
extensions of the Higgs sector. The propagation of quarks of masses 
Q>^ ■ up to about 5 GeV are considered, taking into account their markedly 

(— I , different dispersion relations due to propagation through the hot elec- 

QhI troweak plasma. It is shown that the contribution of these lighter 

O ; quarks to the baryon asymmetry can be comparable to that for the t 

D . quark considered earlier. 



I. INTRODUCTION 

It has been an extremely interesting observation that all the conditions 
needed for baryogenesis [1] are already present, in principle, in the standard 
electroweak theory. The baryon number violation in this theory, although 
exceedingly suppressed at the present time [2], can be unsuppressed at high 
temperature [3]. C and CP violations are contained in the interaction of 
quarks with Higgs fields. Finally departure from thermal equilibrium, al- 
though difficult to obtain at the electroweak scale - typical weak interaction 
rates are extremely faster than the expansion rate of the Universe - can 
nevertheless exist if the electroweak phase transition is of first order. We 
thus have the exciting possibility of explaining one of the most fundamental 
problems of cosmology in terms of laboratory physics. 

Several mechanisms have been proposed to obtain the observed baryon to 
entropy ratio in the electroweak theory [4-9]. In particular. Nelson et al [7,8] 
consider simple extensions of the minimal standard model (MSM) to obtain 
sufficient CP violation. This is provided by the complex space dependent 
fermionic mass function within the bubble wall, which arises quite generally 
in such models. The reflection and transmission coefficients are then different 
for particles and antiparticles, leading to a separation of some CP-odd charge. 
The latter is then converted to baryon asymmetry in the broken phase by 
the baryon number violating process in the unbroken phase. They consider 
only the propagation of the top quark through the medium due to its large 
Yukawa coupling [8]. 

The MSM, whose CP violation was earlier thought to be too small to gen- 
erate any significant baryon asymmetry, has recently been shown by Farrar 
and Shaposhnikov [9] to have the potential to generate this asymmetry by 
the above mechanism, provided one takes the quark mixing effects at high 
temperature properly into account. Further they consider a direct separation 
of the baryon number by the bubble wall rather than of some other CP- odd 
charge. 

Although the nonminimal models were initially studied because of the 
insufficiency of baryon asymmetry produced in MSM, we believe that it is 
worthwhile to investigate these extensions on the cosmological front, as long 
as one is looking for deviations from the MSM in the laboratory [10]. 

Here we reconsider the simply extended MSMs for lighter quarks with 
masses up to about 5 GeV. The dispersion relation for these quarks in the 



electroweak plasma is markedly different from the free one [9,11]. It is not 
immediately clear how such propagations are going to affect the calculation 
of the baryon asymmery. Following Ref [9] we take the altered fermion prop- 
agation into account to leading order and examine the generation of this 
asymmetry by the so-called normal and abnormal modes. 

We avoid the details of specific models by parametrising the mass func- 
tion in a simple way and consider the direct separation of baryon number 
by the wall. We also restrict our calculation to low bubble wall velocity. 
The reflection coefficients of quarks from the bubble wall are obtained in an 
iteration series. We argue that our calculation, though oversimplified, does 
give the order of magnitude estimate of baryon asymmetry, subject to the 
uncertainties involved in the details of phase transition and other properties 
of the medium. 

In sec. II we review the propagation properties of quark excitations in the 
electroweak plasma. We then calculate in sec. Ill the reflection and trans- 
mission currents giving rise to baryon asymmerty. These expressions are 
then evaluated in sec. IV. Our concluding remarks are contained in sec.V. In 
the appendix we solve the quark equation of motion within the wall in an 
iteration series to fourth order. 



II.QUARK PROPAGATION IN HOT PLASMA 

Here we review the propagation properties of light quark excitations in 
the electroweak plasma at high temperature. The most important effect of 
the medium on the quark is that it acquires a chirally invariant effective mass 
with an altered dispersion relation. Neglecting smaller contributions due to 
the weak gauge boson and Higgs scalar exchanges compared to that due to 
gluon exchange, the one loop self-energy leads to the same effective mass Eq, 
for both left (L)- and right (R) -handed quarks, 

Eo = (27ra,/3)^/2T ^ 0.5T (1) 

with 0:5 = .12 at the Z boson mass. For excitations close to Eq, the effective 
Lagrangian incorporating the altered dispersion relation is [9] 

C = iR\do + -(T.V + iEo)R + iL\do--a.V + iEo)L + mL^R + m*R^L (2) 



where we have also included the mass acquired through Higgs mechanism. 

The Lagrangian (2) gives the equation of motion for the L and R com- 
ponents. In the following we consider the one-dimensional problem where 
quarks propagate along the z-axis, normal to the bubble wall. Writing 




the equations split into two independent sets. Defining 





and considering solutions with positive energy E, they are 



where 



^4 = iQ{z)<S>, (3) 



QW=3r .r _(i"_j5j. (4) 




and a similar one for $' with m replaced by its complex conjugate. These 
equations refer to the fluid rest frame. Although we work in the wall rest 
frame, it suffices to evaluate the reflection coefficients in the fluid rest frame, 
because of our restriction to linear terms in the bubble wall velocity in cal- 
culating the baryon asymmetry. 

The planar bubble wall has a finite thickness, extending from 2; = to 
z = zq. It separates the broken phase {z > Zq) from the unbroken phase 
{z < 0). The Higgs induced mass m rises from zero in the unbroken phase 
through the bubble wall to the (almost) zero temperature mass mo in the 
broken phase. 

Requiring plane wave solutions in the broken phase, eqn.(3) yields the 
altered dispersion relation mentioned above [12], 



E^ = Eo±^ip/3y + ml (5) 

In contrast to the free particle dispersion relation, where only one branch 
belongs to positive energy, the presence of Eq{> mo) in (5) now makes both 
branches accessible to a particle with positive energy. The two branches are 



called normal (+) and abnormal (-) modes of propagation. Unlike the energy 
E, the value of the variable p, however, does not give the true momentum 
p± of the excitation, the latter being given by 

Also the group velocities for the two branches are given by 

v± = --E± = ± , ^=. (7) 

dp 3 ^(^73)2 + ml 



In the unbroken phase, where the dispersion relation becomes 

E± = Eo± k/3, (8) 

each of the components ipi satisfies an uncoupled equation. ■?/'i 2 , belonging 
to chirality +1 , has k = ±3{E — Eq) respectively, while ■?/'3^4 , belonging to 
chirality —1, has k = ^3{E — Eq) respectively. It might be thought that 
a particular component ipi would describe a left- or right-moving particle, 
depending on whether the energy belongs to the normal {E > Eq) or the 
abnormal {E < Eq) branch respectively. However, this is not true, as the 
direction of the true momentum (and also of the group velocity ) is opposite 
to that of k in the abnormal mode. 

Over the domain wall, m is space dependent. We may write the solution 
for ^{z) as 

^z) = n{z)^{0), 0<z< zq, (9) 

in terms of the 2x2 unimodular matrix Q{z). At z = Zq its elements will be 
denoted by 

"(^«) = ( ;. f. ) ^ <'•" 

An iterative solution is obtained in the Appendix. 

We note here the Lorentz invariant expression for the density of fermionic 
excitations, 

n = {exp Pp ■ V + 1)^^. 

Here f3 is the inverse temperature of the fluid in the frame where it is at rest, 
pf^ is the energy momentum 4-vector of the excitation and v'^ , the 4-velocity 



of the medium. In the wall rest frame, p'^ = {E,p) , v^ = 7(1, f) where 
7 = l/yl^-t^ and p is the true momentum given by (6). For p along the 
positive z-direction, we thus have in this frame, p ■ v = E±{1 =1= v/3), up to 
linear term in v. In the following we need the densities of particles moving 
towards the wall. In the unbroken phase these are given by 

^± " e/3£±(i-t>/3) ^ X (11) 

for the (±) modes respectively. In the broken phase the corresponding quan- 
tities n\. are given by the same expressions with the reversal of sign before 

V. 



III. REFLECTION AND TRANSMISSION CURRENTS 

Since baryon non-conservation through the sphaleron processes involves 
the left-handed fermions and antifermions, we are interested in calculating 
the left-handed baryonic currents only. 

Consider first the propagation of quark excitation by the normal mode. 
We send a right-handed fermion towards the domain wall from the unbroken 
phase. Noting the revarsal of chirality at the wall, the incident wave (of unit 
current) and the refiected wave of amplitude r , say, is given by 

^z)= (J)e^'=^+(^°)e~^'=^z<0. (12) 

On the right ( broken phase), we have only the transmitted wave of amplitude 
t, say. Solving eqn.(3) for $ we get 

Here p satisfies eqn.(5) with the plus sign and 

cosh^ 



E-E0+ p/3 



\ 2p/3 • 

To find the unknown ampltudes we use the boundary conditions given by 
eqn.(9) for z = Zq, 
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The reflection coefficient is 

+ ' ' |a*cosh^ + /3sinh^|2' ^ ' 

The incident current, same for particles and antiparticles, is |^+, where n^ is 
given by (11). Then the net contribution to the reflected left-handed baryonic 
current is 

/ ^-nl{R^-R^) (16) 

Here and in the following a bar on a reflection or transmission coefficient 
denotes the corresponding quantity for the antiparticle. It is obtained by 
solving the same eqn.(3) with m replaced by m* . 

Next we calculate the transmitted baryonic current in the unbroken phase 
due to incidence on the wall from the broken phase. On the left there is simply 
a transmitted wave of amplitude t, say, 

^z)=(^~\e-'''\z<Q (17) 

On the right we have both the incident wave (of unit current) and the re- 
flected wave of amplitude f, say. Solving eqn.(3) separately for the two cases, 
we get 

^ , , / sinh^ \ ■ / ^ ^ _ / cosh^ \ • /^ , •. , , 

y — cosht^y y — smhfc'y 

Again using the boundary condition (9) at z = zq, we have 



sinh 6* + r cosh 6* \_(a P \ f^ 
-(cosh^ + fsinh^) J ~ \ P'' a* I \i 



(19) 



The transmission coefficient is then 

T+ = |tp = l-R+ (20) 

There arises a transmitted left-handed baryonic current in the unbroken 
phase given by 





We add (15) and (20) to get the total normal mode baryonic current in 
the unbroken phase due to reflection and transmission as 

T r/A-1 
J+= J --{nX-nl){T^-f^) (22) 

3mo 

A similar contribution to the baryonic current results from propagation 
by the abnormal mode in the energy region E < Eq — mo. It is given by 

'r°rfA;l 
J-= I --{nt-nl){T^-f^) (23) 



where 



with 



3mo 



X_ = (24) 

|«*cosh^'-/5sinh^'|2 ^ ^ 



cosh 9' 



Eo-E + p/3 
\ 2p/3 • 



The total CP-violating left-handed baryonic current in the unbroken 
phase generated by reflection from and transmission through the bubble wall 
of the fermionic excitations by the normal and the abnormal modes is 

J^^ = J^ + J^ (25) 

The final step is to obtain the baryonic density hb in the broken phase 
from the steady state solution to the rate equations in the two phases. Nelson 
et al find numerical solution to the Boltzmann equations. We shall follow 
Farrar and Shaposhnikov, who solve the diffusion equations for small bubble 
wall velocity to get 

UB = JSpf (26) 

where / is a given function of the diffusion coefficients for quarks and leptons, 
the wall velocity and the sphaleron induced baryon number violation rate. 
Their estimate for / is 10^^ < / < 1 in MSM, which should also be valid for 
its simple extensions. 



IV. ESTIMATION 

In the standard model with a single Higgs doublet, the expectation value 
of the Higgs field is real everywhere during phase transition. Adding extra 
multiplets will allow in general some of their components to acquire space- 
dependent complex values within the bubble wall. This in turn leads to 
a complex mass function for the quark having Yukawa coupling to those 
multiplets. It is, in principle, derivable from the model considered but, in 
practice, will depend on the (unknown) Higgs couplings. Here we assume the 
simplest form for the mass function, 

m{z) = — z + i—z{zo — z) (27) 

within the bubble wall. The parameter 5 is related to the magnitude of CP 
violation in the model. 

With the parametrization (27), we solve eqn.(3) by iteration. This solu- 
tion becomes a power series in three dimensionless quantities, viz, the energy 
variable y = 3{E — Eq)zo and the constants c = 3mQZo and d = 36zo. In the 
Appendix we have obtained this series up to fourth order. Inspection of the 
coefficients suggests that for y, c and d less than unity, it should represent 
the solution well. 

The condition y < 1 restricts kto k < 1/zq. Now the refiection coefficients 
are known to be asymptotically ~ exp{—27i\k\zo) which appears to set in 
already for k > 1/zq. Our approximation then consists of truncating the 
upper limit of k integration at A; = 1/zq, within which we use our solution to 
compute the reflection coefficients. 

As long as the wall thickness zq < {l5GeV)^^ [13,14], the condition c < 1 
would allow quarks with masses up to that of the b quark. Note that even if 
this condition had allowed the t-quark mass, our result would not apply to 
this quark propagation, since the dispersion relation on which we base our 
work, would not be valid. Instead, the free particle dispersion relation would 
be more appropriate for the t-quark, as has been assumed in the work of 
Nelson et al. Finally the condition d < 1 or 6 < l/3zo is a very reasonable 
bound for the imaginary part of the quark mass within the bubble wall. 

It is now simple to evalulate the currents J± given by eqns. (21,22). In the 
following we choose Zq ~ {20GeV)~^, whence the upper limit of k integration 
is ~ 20GeV. As the temparature of the phase transition ~ lOOGeV^, we may 
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expand the density function in v for small v and approximate the exponentials 
by unity to get 

nl-nlr:^-\(3E^v. (28) 



The transmission coefficient may be calculated in a straightforward way using 
the values of a and jS given in the Appendix. We get 

T -f - Tt-^r v'^^^^(V3 - 4I/V45 + 5cV42) 

^ ^ ^ {(l-cV6-cVl68)i/ + cV/45 + v/^^^^}2 ^ ' 

Clearly for c < 1, the higher powers of c can be safely neglected. Thus 



^^ = ^^7^ L ^^ {y + Vy^ - c^r ^ ^ 

Observe the large cancellation in the sum of J+ and J_ arising from quark 
propagation by the normal and abnormal modes respectively. The asymme- 
try current becomes 

,^ _ 2/3vdc /•! y^/]p^^{l - Ay'^/15) 



ISvr 



l3v5mQ (31) 



Finally noting the one dimensional entropy density s = 737r/3/5, the baryon 
to entropy ratio is given by 

ub/s ~ lO-%f5mo (32) 

where 6 and rriQ are expressed in GeV. Within the framework of our cal- 
culation, it is the dynamics of the b quark propagation through the phase 
transition bubbles, which gives a sizeable contribution to hb/s. With v ~ 0.1 
and 6 not too small, its contribution can be comparable with the observed 
asymmetry, hb/s ~ 5 x 10"^^. 
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V. CONCLUSION 

We have investigated the generation of baryon asymmetry in the standard 
model of the electroweak theory with more than one Higgs multiplet. Such 
models generally give rise to a complex mass function for a quark within 
the wall of bubbles formed during the phase transition. This constitutes the 
CP violation needed for baryogenesis. We do not attempt to calculate such 
a mass function, however: We paramertise its real and imaginary parts in 
a simple way, making the integrals trivial to evaluate. We follow Ref [9] 
to consider the direct separation of baryon number by the phase boundary 
rather than separation of some other CP-odd charge to be converted into 
baryon asymmetry by a separate process, as discussed in Ref [8]. 

The inclusion of the temperature dependent effective mass gives rise to 
two modes of quark propagation in the plasma. Our calculation shows that 
both modes must be considered. In fact, the net baryon asymmetry cur- 
rent results after large cancellation between the baryonic currents carried 
separately by the two modes. 

The calculation presented here is a simple analytic one giving the order 
of magnitude estimate of the baryon asymmetry. Our formula (32) neatly 
isolates the model dependent parameters involved in the description of the 
electroweak medium and the first order phase transition. These are the wall 
thickness zq [15], the velocity v of the medium, 6 giving the magnitude of the 
imaginary part and / related to the plasma diffusion and sphaleron transition 
rate. 

The calculation is limited to small bubble wall velocity in the rest frame 
of the plasma, as we keep only the terms linear in v. Also it concerns lighter 
quarks for which the dispersion relations deviate appreciably from those of 
free propagation. For the t quark the free particle dispersion relation is 
accurate enough and its propagation can be discussed following the treatment 
of Nelson et al. Our calculation indicates that the contribution of lighter 
quarks, like the b quark, to the baryon asymmetry can also be substantial. 

After completing the work, we learnt of several works [16,17] taking into 
account the effect of damping in the quark propagation. These authors have 
shown that it reduces the reflection coefficient to negligible values making 
the present mechanism totally ineffective to reproduce the observed baryon 
asymmetry in the minimal standard model of the electroweak theory. How- 
ever, this conclusion is not agreed by others [18], who claim that the damping 
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rate is irrelevant to the problem. 

Even if we take the damping into account, it is easy to see that it cannot 
drastically affect our result obtained in nonminimal models. In the mini- 
mal model, the baryon asymmetry formula involves high powers of the mass 
matrix. Damping effectively multiplies each such matrix by a suppression 
factor, making the asymmetry negligible. However, in our formula the mass 
function appears only quadratically. Furthermore, the range of values al- 
lowed for the wall thickness zo is expected to be large in nonminimal models. 
For zq < {20GeV)^^, a simple estimate indicates that the total suppression 
factor cannot be smaller than 10^^. Inclusion of this factor still keeps our 
result (32) for hb/s near the observed value, given the large uncertainty in 
the value of / mentioned at the end of sec. III. 

One of us (S.M.) gratefully acknowledges the hospitality at the Institute 
for Theoretical Physics, University of Berne, where this work was started. 
He also thanks Prof. M.E.Shaposhnikov for an useful discussion. 

APPENDIX 

We solve eqn. (3) by iterating the corresponding integral equation, 

^z) = $(0) + I r dz^Q{z^)^z^) (33) 

Jo 

For the matrix Q{z) in eqn. (9) the solution reads 

pZ pZ P Z~i 

Q{z) = l + i dziQizi) + e \ dzi / dz2Q{zi)Q{z2) + ■■■ (34) 

Jo Jo Jo 

It is convenient to express Q{z) in terms of Pauli matrices a"^{m = 1, 2, 3), 

with 

/i = -ilm m, /2 = iRe m, f^ = E - Eq 

With repeated use of 

we can reduce the products of cr™ to a linear combination of o"™" and the unit 
matrix. Writing 

n{z) = UJo{z)l + uJm{,z)a"' 
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we obtain the iterative solution for the matrix elements up to fourth order, 

fZ rzi 



rz rzi 

uio{z) = 1 - / dzi dz2i{zi) ■ f{z2) 
Jo Jo 



+ r dz, r dZ2 r dz,e"'^Pfmiz,)Uz2)fpizs) 

Jo Jo Jo 

rz rZi J-Z2 r-zs, 

+ I dzi I dz2 I dzs / dz4{f{zi) ■ i{z2)i{zs) ■ f (-24) 

- f (zi) ■ f(^3)f (^2) ■ f (^4) + f (^1) ■ i{z,)f{z2) ■ f (^3)} (35) 



rz rz rzi 

Ura{z) = i dZifra{Zl)-i dZi rfZ2e"'^'/fc(^l)/«(^2) 

JO Jo Jo 

rz rzi rZ2 

- i dzi dz2 / dzs{fm{zi)i{z2) ■ f (23) - fm{z2)f{zi) ■ i{zs) 
Jo Jo Jo 

+ fU^s)f{z,) ■ f (^2)} + ^ [ dzi £ dz2 fj dz, £ rf^4[e'"^^ 

{/p(^l)/,(^4)f (^2) ■ f (^3) - /p(^2)/,(^4)f (^1) ■ f (^3) 

+ /,(^3)/,(^4)f (^1) • f (^2)} + e^''fa{z^)Mz2)Uzs)U^,)], (36) 

Here f(-2i) ■ f(-22) = fm{zi)fm(yZ2), for example. With the parametrization 
(26) it is easy to evaluate the integrals. We actually need a;o,m(-2o) which we 
simply denote by a;o,m-Using the variables y, c and d introduced in the text 
and rejecting powers of 6 higher than the first, eqns. (32,33) give the following 
expressions {u = iu^) 

ujo = 1 y^ -\ \ y* y^ -\ 

" 2^ 8 24^ 720^ 

+di-^y + ---) (37) 

c Co 3c^ 

c c 2 c^ 

^(^,y^-) (39) 
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The matrix elements of ^{zq) in eqn.(lO) are given by 

a = Wo + iuJs, 13 = uoi- iuo2 
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